Spherical spin-glass— Coulomb gas duality: solution beyond mean-field theory 
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We present an alternate solution of a Gaussian spin-glass model with infinite ranged interactions 
and a global spherical constraint at zero magnetic field. The replicated spin-glass Hamiltonian is 
mapped onto a Coulomb gas of logarithmically interacting particles confined by a logarithmic single 
particle potential. The precise free energy is obtained by analyzing the Painleve r^^[n] function 
in the n ^ limit. The large A'' thermodynamics exactly recovers that of Kosterlitz, Thouless 
and Jones [J]. It is hoped that the approach here can be extended to apply to systems beyond the 
spherical model, particularly those in which destabilizing terms lead to replica symmetry breaking. 
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Spin-glass systems are an important component of the 
longstanding effort to understand the statistical mechan- 
ics of disordered systems. In particular, spin-glass models 
exemplify the interplay of frustration — an outgrowth of 
quenched randomness — and thermal fluctuations ^ . At- 
tention has centered on the phase transition from a high 
temperature paramagnetic phase to the spin-glass phase, 
in which the spins are frozen in random directions over 
macroscopic time scales 0-0] • 

Several decades ago, an alternative yet simple, non- 
Ising spin-glass model was proposed by Kosterlitz, Thou- 
less and Jones (KTJ) based on the exactly solvable spher- 
ical model (SM) [jy] , which has been termed the spherical 
spin-glass model (SSGM). KTJ constructed a solution 
based on the known methodology of the SM that depends 
on the density of eigenvalues of the exchange coupling 
matrix For the SSGM, the exchange coupling matrix 
in the large N limit is effectively a Wigner- Dyson random 
matrix that follows a semi-circle law distribution. One 
of the most appealing features of the KTJ analysis was 
its prediction of a spin-glass transition without the need 
to address replica symmetry breaking, as the quenched 
averaged free energy can be computed without requiring 
replicas at all. A second approach attempted to work 
within the replica framework, pointing out a Parisi like 
order parameter and possible replica symmetry breaking 
perturbations [ol [lo|. 

In this letter, we introduce a third approach that solves 
the SSGM without recourse to the large N random ma- 
trix analysis of KTJ. The SSGM is mapped onto an ef- 
fective Hamiltonian of logarithmically interacting charges 
in a logarithmic single particle potential, the number of 
replicas corresponding to the number of particles n. This 
approach does not require any assumptions about the 
replica structure, and the n — > limit can be extrapo- 
lated using inductive diagrammatic arguments to connect 
the grand canonical and canonical ensembles. The par- 
tition function is shown to be equivalent to the Painleve 
[n] function, which is solved explicitly to yield a spin- 
glass transition that is in complete agreement with the 



results of KTJ. It is hoped that the spin- glass-Coulomb 
gas connection and the methods described here are gen- 
eral enough to be applied to solve both the Sherrington- 
Kirkpatrick and Edwards- Anderson models beyond the 
limitations imposed by existing mean-field and numeri- 
cal approaches. We note that a general relationship be- 
tween replica based Hamiltonians and Painleve transce- 
dents has been rigorously demonstrated ll[ and there- 



fore the approach can be possibly extended to the broader 
category of replica based or disorder-induced critical phe- 
nomena. 

The SSGM is based on a strictly Gaussian model, in 
which the spins take on continuous values in a specified 
range. The Hamiltonian for N spins is 



J ij Si Sj 



(1) 



where the quantity A is the Lagrange multiplier that en- 
forces the spherical constraint, i^f) — ^ ^^'^ each 
exchange coupling J^- is subject to a Gaussian distribu- 
tion, P(Jjj) cx e-^4V2>/'^ ^j^j^ J ^ j^j,^ gj^j,g 
interested in the quenched average of the free energy, we 
apply the replica technique and consider the annealed av- 
erage of the n"^ power of the partition function Z^. In the 
first step we integrate over the Jij 's in the usual manner 
to yield an effective four-spin interaction that couples dif- 
ferent replicas, with a, b replica indices that range from 
to n. Next, these terms can be decoupled by intro- 

ducing an auxiliary field Q, in a Hubbard-Stratonovich 
transformation. 



Zn = exp 



We can then integrate out the S'^'s: 
Zn = exp 



Trln( / A+Q 



N ^2 ■ 



(2) 
(3) 



In the final manipulation, we have to integrate over the 
matrix Q, which is real and symmetric. We are led to the 
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FIG. 1. 1/5, a particular 5-vertex diagram in the virial expan- 
sion of Zn 



evaluation of a partition function over the eigenvalues of 

Q, \i having the form, Z„ = / exp \H.cg{}^\ dXi . . .dA„, 
where the Hamiltonian Hcg corresponds to a one- 
dimensional gas of particles with logarithmic interac- 
tions, which we term a Coulomb gas (CG) [l^ . 



with the single particle potential 



y(A) = ^ln[A 



(4) 



(5) 



For the rest of this paper, we fix A = 2 at the unitary 
ensemble, as this allows us to adapt results now in the 
literature 12- We believe that this substitution pre- 
serves the essential features of the model. 

As a prelude to the full calculation of the partition 
function, we outline the development of a diagrammatic 
approach to its evaluation. In the replica formalism, the 
limit n — >■ must be taken with care. Our basic strategy 
is to first work within the canonical ensemble and observe 
the general the trend of its n dependence in order to 
transform to the grand canonical ensemble, creating an 
adjustable n that can be extrapolated to zero. A useful 
method for treating the logarithmic interaction is a virial 
expansion of 



^ I dXid\2 ■ ■ ■ dXnC 
= ^/dAirfA2---dA„e-^[^ 



-l^JdX^dX^... dA„e-^[^'l (- A In I A, - A, 



Take for instance an arbitrary five vertex diagram, that 
we term L5 as shown in Fig. [1] The expression to which 
this diagram corresponds to is given by: 



L5 = M5 J dXidXjdXkdXidX„^e ' (- 
e (— Aln |Afc — Xj \) e x 
e (-Aln|Am - Ai|)e 



■Aln|Aj - Ail) X 



(7) 



where the quantity = n\nin2/{{n — 5)!3!2!) is the 
overall weight of the diagram. Next, one accounts for 
the number of ways of assigning the n total vertices to 
the 5 vertices in the diagram. Then, one multiplies by 
symmetry factors from permuting lines connected to a 
common vertex, yielding ni = 712 = 1/2. One straight- 
forwardly infers that any diagram is simply a product 
of expressions involving connected diagrams followed by 



the factor Z, 



(0) 

n-M 1 



where the "bare" partition function is 



given by, Z 



(0) 



jdXe 



-V[X\ 



The grand partition function, 0(z) — -n -^n^^ can 
be formally constructed with each vertex multiplied by 
a fugacity factor z such that Q'^'^\z) — exp[z(7o], where 
go = J dXe . In the limit n 00, Q{z) and Qn 
satisfy the relations, lnZ„ = lnO(z(n)) — nz{n) and n = 
^8(z). To extract the rt — > limit, the n dependence 
of L5 reduces to n x 4!(— 1)''. Dividing out n, and re- 
expressing the factorial in terms of a F function, we are 
left with the remainder — (— 1)^ /q t'^e~*dt. 

Next we make use of the zeroth order, non-interacting 
term in the partition function, « 1 -|- nlngo- Since, 
L5 scales as 1 /q^ , we can absorb this factor by recasting 
L5 as, -{-If dt{t^e-i°*/t). Finally, this allows the 
virial expansion to be completely factored in terms of 
its n dependence. The integral form of the logarithm. 
In (70 = Jq rfi(e~* — e~**)/i can be substituted with 



Q^^^t), or the non-interacting grand partition 
function at negative fugacity. The disorder-averaged free 
energy is, then. 



(InZ) 



Additionally, if the 
n-component system 

easily show that Eq, 
replica calculation 



dt{e~^ - e{-t))/t 



(8) 



partition function of the 



takes 



the form, Z„ = 
• ■ ^k) Y[i=i dxi then one can 
© is equivalent to the standard 



(InZ) 



lim (Z„ 



1) /n = dZn/dn\^ 



(9) 



Subsequently, the equation for the density can be repre- 
sented as a diagrammatic series where p[x) = poX K and 
the graphical form of K is shown in Fig. [21 K can be 
generated by observing that V{X) couples to the density 



(6) p at the zeroth order, po(A) 



ze 



The tree-level 



diagrams of Fig. [2] (a) result in a self-consistent series 
for mean field density pmf{x). The summation of these 
terms reduce to an inhomogeneous Sine-Gordon equa- 
tion. The proper solution is non-trivial and a thorough 
analysis will be discussed in a separate article. In princi- 
ple one should recover the known KTJ result: 
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T<Tc 
T>T, 



(10) 
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(a) 0+0- 



(b) 



K= (a) + (b) 

FIG. 2. The full diagrammatic expansion of K. An empty 
circle stands for the vertex at point x and the dark circles are 
at other points. A line is simply the distance function \x — x'\. 
(a) The tree-level diagrams resulting in the mean-field density 
Pmf{x) (b) Loop contributions beyond mean- field theory. 



We now leave our discussion of the mean field theory 
and focus on the full solution, which can be evaluated 
exactly using r function theory of Painleve systems. The 
first task is to recast the Okamoto T''^^[ n] 14| integral 
discussed by Forrester and Witte (F&W) [13[ into a form 
that is equivalent to the partition function of Eq. (|4]), 



dxi--- 1 dxn Jl e~^3 {t - XjY W {xk ~ Xj) 

-oo J-oo l<J<k<n 



(11) 



The integrand above can be rescaled by letting Xj 
yj\/Na, and let t —J- A\/Na, yielding 



Z„ = (7Va)("(«-l)/2+n+M)/2x 



l<j<k<n 



Vj) 



(12) 



and the n — 1 term is given by, 

t'^[1] = / [t-xYe-'^'dx 



(13) 



Recurrence relations for t^^ [n] have been thoroughly in- 
vestigated by F&W in the context of certain random ma- 
trix averages. These relations are well suited for our pur- 
poses, 



a^^[n] = det 



r7^^[l] 



(14) 



where the determinant is of an n x n matrix with indices 
(i,j) and the relationship between and r ''^[n] is 

given by. 



' ri-l 



= 2"("'i)^"/2 I n I "r^^H (15) 



It follows that the n = 1 part 



jv 



1] = 



1] 



(16) 



is used to generate higher orders in n via Ea. (jl4p . This 
allows us to extrapolate the n dependence of T^^[n] in 
the case in which /i = —N/2 for large postive N . Let's 
place t slightly off the real axis and consider the integral 
of Eq. ([T^ when /i = — 1, which can be shown to take 
the general form: 



{x-ty\^'\e"'^dx 



-l)!d<lA'l-i V Jo J 



(17) 



which is equivalent to the Dawson integral 15 1 



dk 



sin(H) dk (18) 



The evaluation of the derivatives of the integral as given 
by the R.H.S. of (fT8)) can be now taken. 



dt\^'\- 



e sin(fct) dk 



— Im 
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/ il'^l-i exp [-k^ + (1^1 - 1) In fc ikt] dk 
Jo 



We can apply the stationary phase approximation and 
evaluate the integral by expanding about the cxtrcmum 
to yield the following result. 



Vt^-8(|p|-l) 



T'-[l]= 

exp 



8 



:H-l)ln((i^^)-^-| 

(19)' 

Equation ()19p can be simplified by the functions f{q) and 
g{q), with q = t/ \/\^\ — 1 ~ t/\/N such that. 



T^^l] =5(9)cxp[7V/(<z)] 



(20) 



We find empirically that the precise leading order in n 
contribution to the determinant in (|14p yields, 



^^'''[n] = I n I .9(9)"e^"/^(')/r(g)"("-i)/^ (21) 

where fi{q) = f{q) + q^ . The final expression becomes, 

r^^[n] = e-*'"2-"("-i)^-"/25(g)"e"^-^'(«)/"(g)"^""'^^' 

(22) 

Taking the appropriate derivative of Eq. ^ yields the 
n — )■ limit. The prefactor of Eq. (IT^ contains no inter- 
esting thermodynamic information; the key contributions 



4 



to the free energy become 
J-[A] = 

-r(lnZ[A])=T(iV/(y) 



In 



2.g(y) 



(23) 

The Lagrange multiplier can be eliminated by enforc- 
ing the spherical constraint via the relation — ^^'"3^^^^^ = 
Ef=i (5?) = N, yielding: 

i(iVr/J) (- Vr2A2/j2 _ 8 + TA/jj = N (24) 

having a solution of A> = 1 + where the range of T 
lies above the branchpoint structure of the squareroot. 
However, for low enough values of T, there is no proper 
mathematical solution, and A< — ^^^+0{1/N) in order 
to obey the spherical constraint. Hence, the complete free 
energy is given by. 



^=-r(lnZ[A>]) 

NT ( (2+T'") 



+8 In 



_2_ 

T' 



AT' 



-4 



(25) 

where T' = T/J. The critical temperature is obtained 
from the maximum of Eq. p5|) , yielding Tc = \f2J . After 
taking the large N limit, we can compare the results to 
KTJ. Taking the appropriate derivatives of T we have 
for the specific heat per site. 



Cv{T) = 



1 



2 - {T.lTf 



T <Tc 

T>Tr. 



(26) 



in precise agreement with KTJ, as the phase transition 
contains a discontinuous 3rd derivative or specific heat 
exponent a = —1 similar to that observed in the three 
dimensional ideal Bose gas and the SM. Although our 
analysis is taken in the unitary ensemble (A=2) rather 
than the orthogonal ensemble (A=l) appropriate to the 
model considered by KTJ, we expect no great discrep- 
ancies as we can generalize our results to account for 
the possibility of other ensembles by including an ad- 
ditional factor B in the spherical constraint equation. 
Re-expressing Eq. 



BN = 2A^T'^ A - VA2 - T'-^ 



(27) 



and solving it explicitly, we find that it is satisfied if T' > 
B/2 and therefore we expect that B and subsequently 
the ensemble characteristics should only affect the value 
of T,. 

To conclude, we should emphasize certain key points 
regarding importance of the mapping utilized here and its 
application to other spin glass models. First the mapping 
follows from a reorganization of the total phase space 



including the replica degrees of freedom without neces- 
sitating any ansatz of the Q matrix structure (replica 
symmetric or otherwise) and thus only the ensemble of 
the disorder is chosen. Subsequently, within the CG rep- 
resentation, a direct confirmation of a thermodynamic 
phase transition can be observed without any knowledge 
of the proper order parameter (droplet or Parisi-like) of 
the original spin glass model. Furthermore, a virial ex- 
pansion of the CG density allows in principle for the sepa- 
ration of mean-field (tree level diagrams) and correlation 
effects (loop diagrams) that may lead to replica symme- 
try breaking. Within the CG system of the short-ranged 
Ising spin glass systems in finite dimensions, the mean 
field theory should precisely recover the Sherrington- 
Kirkpatrick solution, and the correlation effects can be 
treated as perturbations which can affect the stability of 
the spin glass phase. 

In the case of the SSGM, the facility in deriving the 
the CG model results from the integrability of the spin 
variables, however for Ising-like spins the methods pre- 
sented here must be replaced by a different manipulation. 
Nevertheless, our emphasis is placed on the known rig- 
orous mathematical understanding of replica based mod- 
els, for which the trace over the total phase space can be 
generically re-expressed as some type of a CG ensemble 
properly taken in the n — > limitjllj. Thus, a success- 
ful application of this approach has been achieved here, 
and we look forward to developing the methods further 
in other spin glass models. 
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